Coherent Spin-Current Oscillations in Transverse Magnetic Fields 
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We address the coherence of the dynamics of spin-currents with components transverse to an 
external magnetic field for the spin-1/2 Heisenberg chain. We study current autocorrelations at 
finite temperatures and the real-time dynamics of currents at zero temperature. Besides a coherent 
Larmor oscillation, we find an additional collective oscillation at higher frequencies, emerging as a 
coherent many-magnon effect at low temperatures. Using numerical and analytical methods, we 
analyze the oscillation frequency and decay time of this coherent current-mode versus temperature 
and magnetic field. 
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FIG. 1. Quasi-classical sketch of transverse spin transport in 
a spin chain directed along the y-direction: (a) External mag- 
netic field with a static 'bias' component B z and a perturb- 
ing, space- and time-dependent component B x (y,i); (b) B" 
produces a magnetization M with a homogeneous transverse 
component and a Larmor precession, B x {y,t) produces an in- 
homogeneous transverse component driving (c) a transverse 
spin-current j x (y,t) in y-dfrection. The dynamics of j x (y,t) 
is controlled by the intrinsic exchange coupling J. (d) As a 
function of time, the magnetization relaxes and the transverse 
component dephases. 



Controlling quantum coherence is paramount for fu- 
ture information processing [l| . The coherence of local- 
ized quantum spin degrees of freedom has been studied in 
a wide variety of systems, including semiconductor quan- 
tum dots [H, Q , molecular magnets Q .nitrogen vacan- 
cies in diamond Jj|, carbon nanotubes H, and ultracold 
atoms Q. Coherence in spin transport has been ad- 
dressed primarily in semiconductors [8j . A new route into 
coherent spin transport may arise from quantum mag- 
nets. Here, magnetization is transported solely by virtue 
of exchange interactions and (de) coherence of spins will 
emerge as a purely intrinsic many-body phenomenon. In 
one-dimensional (ID) spin systems, magnetic transport 
has experienced an upsurge of interest in the last decade 
due to the discovery of very large magnetic heat conduc- 
tion Q and long nuclear magnetic relaxation times fioj ]. 
Genuine spin transport in such materials remains yet to 



be observed experimentally and if combined with mate- 
rials with small exchange couplings [ll], [l2[ , the coherent 
manipulation of spin transport using laboratory magnetic 
fields may become feasible. Theoretically, spin transport 
has already been given significant attention (see Refs. [l3| 
andllj for a review), previous studies, however, have fo- 
cused on the longitudinal spin conductivities only, ex- 
cluding the physics of coherence. Therefore, in this Letter 
we investigate the dynamics of spin-currents with compo- 
nents transverse to an externally applied magnetic field, 
as sketched in Fig. [T] This setup allows us to study 
the collective precession frequency of the transverse spin- 
current and its decay time, which will be at the prime fo- 
cus of this Letter. We will show that, besides a coherent 
oscillation at the Larmor frequency, a second nontriv- 
ial collective oscillation at higher frequencies emerges at 
low temperatures. This oscillation is identified as a pure 
many-magnon effect and also becomes coherent in the 
low temperature limit. 

In this Letter, we study the antiferromagnetic Heisen- 
berg spin chain, one of the fundamental models to de- 
scribe magnetic properties of interacting electrons. It is 
relevant to the physics of low-dimensional quantum mag- 
llal. ultra-cold atoms [Iff, nanostructures 171, and 
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seemingly unrelated - fields such as string theory [18 1 



and quantum Hall systems [19| . The Hamiltonian reads 
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5° (a = x, y, z) are the components of spin-1/2 operators 
at site r, N is the number of sites, J > is the exchange 
coupling constant, B z is a longitudinal magnetic field, 
and H = 1 For B z < Bl = 2 J, Eq. (TJ implies 

a gapless Luttinger liquid [2l|, [22j and, for B z > B*, a 
gapped ferromagnetic ground state. 

We investigate the transverse spin-current dynamics 
for two complementary scenarios and use methods appro- 
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priate for each situation. First, we study current autocor- 
relations at finite temperature, using numerically exact 
diagonalization (ED) and an asymptotic analytic analysis 
(AAA). Second, applying the adaptive time-dependent 
density matrix renormalization group (tDMRG) (23[, we 
analyze the real-time dynamics of currents at zero tem- 
perature during the evolution from initial states with an 
inhomogeneous magnetization. Qualitatively, the same 
physics explains our observations in both scenarios, and 
even a quantitative agreement can be obtained. 
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P E ™{l\j»\m){m\j v \l)6(u>-E m + 
\m) and Ei >m are eigenstates and 
: 1/T is the inverse temperature, 



We begin by discussing the current autocorrela- 
tions C^ja;)=£ lTO e 
Ei)/ZN [241, wh ere \l) 
-energies of Eq. ([T]), j3 

is the zero-momentum spin current, /i, v = x, y, z, Z 
is the partition function, and uj is the frequency. More 
precisely, we consider a symmetrized version C^ v (uj) = 
[C^ v (oj) + C fi "(~uj)}/2, i.e., in the time domain, we fo- 
cus on the real part C» u {t) = Re[C r '"'(t)]. While B z 
breaks total spin conservation, a spin-current can still 
be defined by decomposing the time derivative of the 
spin-density at momentum q as d t S^ = d t S£\j = o — tqj£ 
into a local source (sink) term due to B z (present with- 
out any exchange interactions) and the actual exchange 
mediated spin-current jg. The latter then derives from 
the continuity equation for at B 2 = 0. In turn 

j = iJ2 r Sr x where j = j q= Q. The eigenstates 

(energies) are classified according to the total spin z- 
component J^S^IO = M\l). Since \l) is independent of 
B z and is diagonal at B z = 0, it will remain diagonal 
at any B z . Moreover, by symmetry C xx = C vv . How- 
ever, since j x mediates transitions between sectors with 
AM = ±1 while j z conserves M, the autocorrelations 
C xx (transverse) and C zz (longitudinal) differ at B z ^ 0. 
This difference is solely due to the field dependence of the 
eigenenergies. Formally speaking, this aspect is at the 
center of this Letter. For the remainder we abbreviate 
(jxx(zz) ky (jx(z) _ ]\j ^ e ^hat \yy ^} ie continuity equation 

the current autocorrelations at small q are related to the 
dynamic spin structure factor S^ v (q, uj), exhibiting a ma- 
trix symmetry identical to C^(uj) at B z ^ 0. 

Generically, the longitudinal autocorrelation decom- 
poses into a Drude weight D z and a regular part, i.e., 
C z (uj) = D z S(uj) + C z cg (uj). A significant body of ev- 
idence in favor of D z (f > 0) ^ for B z < B z has 
been gathered EH, 0] , with thermally activated be- 
havior of D Z (T) for B z > B z Ql- Less is known 



on the specific shape of the regular part [26[. The 
previous discussion of symmetries of C^ v implies that 
C x (uj) = T,±[ DXS (u ± BZ ) + C x cg (uj ± B z )]. In general, 
D x and C x cg will not be identical to D z /2 and C z eg /2, 
respectively, due to the different independence of Boltz- 
mann weights in C z and C x . In the time domain, D x im- 
plies a coherent, nondecaying oscillation of the transverse 
current at the Larmor frequency wj, = B z , permitted by 
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FIG. 2. (color online) Frequency and time dependence of the 
autocorrelation C x for B z / J = 2 and (a), (b) f3 = 0, (c), (d) 
j3J = 2 (ED). The Drude weight in (a), (c) is visible exactly 
at the Larmor frequency oil = B z (vertical dashed-dotted 
lines). In (b), (d) the envelopes of fits [as defined in Eq. ©] 
to data for N = 18 are shown (dashed-dotted curves). 



the integrability of Eq. ([T]) [l3( . On the other hand, C x eg 
a priory implies only decoherence and damping. In the 
following, however, we demonstrate that at sufficiently 
low T and finite B z , out of C x , a new collective quasi- 
coherent oscillation of the current emerges. The oscilla- 
tion frequency differs from the 'simple' Larmor frequency 
and cannot be understood within a onc-magnon picture. 

First, we discuss high temperatures, i.e., /? = 0. A 
straightforward analysis yields 



C x {lo) = [C z {uj-B z ) + C z {lo + B z )}/2. 



(2) 



Figure [5] (a) displays C x (oj) for B z / J = 2. As can be 
seen from Fig. [5] (b), this approximately transforms into 
C x (t) ss [R(t) + D x ] cos(i*jLt) in the time domain, with 
R(t) rapidly decaying within ~ 1/ujl and a 'trivial' co- 
herent oscillation due to the Drude weight. 

Next we reduce the temperature to j3J = 2 at B z /J = 
2. Figure [2] (c) clearly shows two effects. First, the 
Drude weight D x is strongly reduced. This reduction 
continues monotonously with increasing B z (as discussed 
within the AAA below) . Second, the regular part C x eg is 
strongly enhanced and undergoes an asymmetric weight 
redistribution with a major peak developing at a fre- 
quency uj larger than ujl (and a minor peak at uj < ujl)- 
This is consistent with B z breaking particle-hole symme- 
try. In the time domain, see Fig. (2](d), we find that 

C x (t) k R(t) cos[(w L + Su)t] + D x cos(uj L t) (3) 

allows for a reasonable leading-order fit of C x (t) over 
several oscillation periods by assuming an exponential 
behavior R(t) = Rq exp(— t/r), i.e., a single decay time r. 
In fact, C x (uj) is rather close to a Lorentzian in Fig.[2](c). 
We find this approach to apply at least to B z / J < 3 and 
(3J < 3 and to have very little finite-size effects for the 
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FIG. 3. (color online) (a) Decay time r and (b) frequency 
shift Suj w.r.t. p for different B z . Data are extracted from 
the autocorrelation C x (t), using ED for N — 18 (symbols). 
Solid lines represent the low-temperature asymptote above 
the critical field, dashed curves are guides to the eye. 



system sizes (N = 10, ... , 18) and time scales (t J < 20) 
studied [2(J. Figure[3](c) is a central result of this Letter. 
It unveils the emergence of a new collective frequency 
scale, namely at ujl + Suj, in the transverse transport 
process besides the Larmor frequency. Moreover, for r — > 

00 this process would be coherent. 

In Fig. [3] we summarize our findings for r and Sui over 
a range of temperatures and fields of < j3J < 3 and 

1 < B z /J < 3. Figure [3] (a) shows that r increases 
roughly linearly with /?, with an increasing slope as B z 
increases. While finite system studies will not clarify if 
this result implies true coherence for a particular range 
of B z as T — > 0, Fig. [3] (a) is at least strongly indicative 
of a large r in that limit. Regarding <5w, Fig.[3](b) clearly 
signals a saturation roughly at f3 J ~ B z / J with the value 
at saturation increasing with B z . We emphasize that 
Suj 7^ directly implies that the transverse current cannot 
be described in terms of transitions between the zcro- 
and one-magnon sector. The dominant spectral weight 
of such one-magnon excitations at q = is exactly at u>l 
25| . leading to Sui = 0. 

To gain insight into the origin of Suj ^ we present an 
AAA for B z > B z and low T. Here, the contribution of 
different M-sectors to C x can be dissected asymptotically 
and, after an extensive analysis [20j , we arrive at a simple 
picture: the transverse current is carried dominantly by 
transitions from the one-magnon sector around q ~ 7r 
into antibound states of the two-magnon continuum at 
the same q. The related frequencies ~ + J and also the 
asymptotic form of C x (t) can be obtained analytically: 
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(4) 



This is consistent with Fig.[3](b), which also suggests that 
Suj — > J as B z increases and (3 — > oo. The thermal acti- 
vation in Eq. ((4]) stems from the one-magnon energy at 
q = 7r. The damping results from summing over all tran- 
sitions in the vicinity of q = tt and its power-law behavior, 
i.e. ~ i -1 / 2 , for t — > oo clearly shows that the single-scale 
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FIG. 4. (color online) Simulation of the transverse spin- 
current dynamics using tDMRG: (a) Initial magnetization 
profiles (S^) at t = for different k; (b) Coherent oscilla- 
tion of the spatially averaged current (jjvyaW) f° r ^ = 0-5 
and B z /J = 1.5; (c), (d) Discrete Fourier transform for (b) 
and for B z / J = 0.75. Besides the dominant peak at the Lar- 
mor frequency Ul = B z , there is another significant peak at 
a higher frequency ujl + Suj. For more details, see Ref. [2(J 



exponential used for R(t) in Eq. ([3]) is an approximation 
for not too low temperatures T only. Nevertheless, for 
a comparison with Fig. [3] (a), we extract a 'decay time' 
from the envelope of Eq. ©, i.e. |Cf sy (i)/C*f sy (0)| < l/e 
for t > t, leading to r = 2/3\/e 4 — 1 ~ 15/3. As shown in 
Fig. |3] (a) (straight line), our ED data for B z / J = 3 are 



consistent with the asymptotic line, e.g., the slope dr/d/3 
is already close to 15/3 at /3J ~ 2. 

Now we turn to the real-time evolution of the spin- 
current derived from a Krylov-spacc based tDMRG ap- 
proach [23|. This allows us to study larger systems than 
with ED, however, at zero temperature and only below 
the saturation field. The latter follows since there are no 
current-carrying states for B z > B z at T — 0. More- 
over, since ED and AAA already suggest that r — s- oo 
as T — >• 0, limitations in the accessible simulation times 
confine the tDMRG to an analysis of Suj. To induce a 
current we add a perturbation Hi = J2 r BfS* to Eq. (P) 
with B'f. = B x cos(27rfcr/7V). First, the ground state of 
H + H\ is evaluated using DMRG, then the system is left 
to evolve under H alone. 

Typical transverse magnetization profiles (S x ) at t = 
are shown in Fig.|4](a) for (B Z ,B X )/J = (1.5, 1) and for 
small values of k. (S x ) follows B x qualitatively, with ad- 
ditional 2fcp-oscillations from the underlying Luttingcr 
liquid. We perform the time evolution using m = 200 
states for the ground-state calculation, a time step of 
St J = 0.25, and a fixed discarded weight [2(J 23 1. Al- 
though Hi breaks £7(1) symmetry, we can still obtain 
reliable results for L < 64 lattice sites. To analyze cur- 
rents free of spatial oscillations we coarse-grain the data 
by averaging over suitable parts of the chain. Figure 
2] (b) shows an example of the time evolution of the cur- 
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FIG. 5. (color online) Frequency shift 5ui w.r.t. the magnetic 
field B\ as obtained from ED (JjJ = 3) and tDMRG. 



rent (i^r/ 2 (*))' averaged over the left half of the chain, at 
(B Z ,B X )/J = (1.5, 1) for k = 0.5. For the times reached 
in the simulation (tJ sa 80), no relaxation can be ob- 
served. This is consistent with the expectation for T = 
drawn from ED and AAA. In view of the substantial 
number of periods covered in Fig. |4] (b) we have chosen 
to directly study the discrete Fourier transform of the 
real-time data in order to obtain the dominant frequen- 
cies. This is shown in Figs. @] (c) and (d) for (J%/ 2 {t)} 
at (B Z ,B X )/J = (1.5,1) and (B Z ,B X )/J = (0.75,1) at 
k = 0.5. These two figures clearly evidence the main 
result from tDMRG, namely that, fully consistent with 
the findings from ED and AAA, there are two charac- 
teristic frequencies in the current dynamics, namely uil 
and ujl + Slo. In contrast to the linear- response regime, 
the analysis of the relaxation scenario finds the contribu- 
tion at ujl to be the larger one. This is not surprising 
since two different scenarios are compared, characterized 
by similar yet different correlation functions. We have 
checked that the results of Fig. 2] (b)-(d) are insensitive 
(i) to the details of the coarse-graining, (ii) to varying k 
within the small fc-regime, and (iii) to the strength of B x , 
at least up to B x /J = 1, as used here [2(J. The latter 
implies a minor role of non-linearity (non-equilibrium) . 

Finally, in Fig.[5]we compare Slo from ED with tDMRG 
for 0.5 < B z / J < 2. The agreement is remarkably good 
indeed, not only in view of the different scenarios. For 
intermediate fields (B z / J = 0.5, 0.75, 1) the frequency 
shifts match each other almost exactly, while we attribute 
the slight deviation of ED from tDMRG at larger fields 
to finite temperature effects, where convergence to the 
T = values has not been reached yet [see B z /J = 3 in 
Fig. [3] (b)]. For B z /J < 0.5, the accessible time scales 
prevent a reliable determination of Slo in our approaches. 

In summary, we studied the transverse spin-current 
dynamics in the spin-1/2 Heisenberg chain. As a main 
result, besides a coherent oscillation at the Larmor fre- 
quency, we provided evidence for a second nontrivial col- 
lective oscillation at higher frequencies, emerging at low 
temperatures as a genuine many-magnon effect and turn- 
ing coherent as the temperature is lowered. 
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SUPPLEMENTAL MATERIAL 

I. Derivation of the Asymptotic Approximation 

In this section we give a detailed derivation of the 
asymptotic approximation in Eq. (4) of the Letter. For 
convenience, we shift the zero point of the energy E to 
N(J/4 - B z /2) and of the quantum number M to N/2. 
The shifted quantities will be denoted by £ and M. in the 
following, 



£=E-N{ 



J B z 



M = -(M-^). (5) 
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Above the critical field the ground state is fully polar 
zed 



i>(M=o,n=o) = ltt •■• tt> 



(6) 



with momentum q = and energy £(o,o) = 0, see Fig. [5] 
(triangle). Periodic boundary conditions are assumed in 
this section. This state is an eigenstate of the current 
operator with the eigenvalue zero 



f V»(o,o) = . 



(7) 



In the limit of T = only matrix elements from this ini- 
tial state contribute to the current autocorrelation func- 
tion C x (lo). Therefore lim^^o C x (lo) = 0. At any finite 
temperature, /3J 3> 1, initial states from the subspace 
with M = 1, i.e. one-magnon states, start to contribute 
to the current autocorrelation function. Using the opera- 
tor T M , which translates a state by \i sites, the eigenstates 
of the Hamiltonian H in Eq. (1) of the Letter read in this 
subspace 

1 N ^ 

q = 2irk/N, k = 0, 1, . . . , N — 1. Their energies are 



£(i,g) = J{cosq- 1) + B z 



(9) 



sec Fig. |6] (crosses) . At non-zero momentum these one- 
magnon states obviously carry a net-current, however, 
they are not eigenstates of the current operator. In fact 
one readily obtains 



J 



f^(l,q) = * IT (1 ~ C tq ) <P(2,q) , 



(10) 



where (p(2,q) are states from the subspace with M. = 2, 
i.e. the two-magnon subspace. They read 
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FIG. 6. The spectrum of the Hamiltonian in the M — 0-, 
1-, and 2-magnon subspaces (symbols). Additionally, those 
transitions of the transverse current are indicated that are 
relevant for the asymptotic approximation at low tempera- 
tures above the critical magnetic field: 'Forbidden' transitions 
(red, dashed lines) yielding no contribution as well as 'allowed' 
transitions (green, solid lines) leading to the dominant con- 
tribution. The latter transitions involve a two-magnon band 
at the edge of the Brillouin zone, i.e., q = n. The associated 
frequencies are larger than the Larmor frequency ujl = B~ , 
namely, these frequencies are shifted by 8u = J. 



In general, and in contrast to Eqs. (JSJ), (HJ, these states 
are no eigenstates of the Hamiltonian H in Eq. (1) of the 
Letter, but 

H v{2, q) = (2B Z - J) <p {q . 2) + i (1 + e««) £ (2 , ?) , (12) 

where <fi(2,q) refers to states generated by the trans- 
verse fluctuations of H , which separate the two adjacent 
flipped spins by one site. However, in the vicinity of 
q = 7r, the prefactor (1 + e lq ) suppresses these contribu- 
tions rendering f{2,it) an exact eigenstate of H. From 
its eigenenergy (2B Z — J) and Fig. [5] (green circle), it is 
obvious that this state is the (anti)bound two-magnon 
state i/j( 2 ^ at q = n, well known from Bethe-Ansatz. Its 



dispersion over the complete Brillouin zone is [27 



S {2,q) = \ (cosg-1) 



J 



2B Z 



(13) 



and is situated above the two-magnon continuum. 

Due to their Boltzmann weight the transitions from 
ip(i.q) into <P(2,n) at g « 7r dominate the current autocor- 
relation function asymptotically for j3J 3> 1. Therefore, 
the leading contribution to C x results from projecting 
all intermediate states solely onto ip r ^ 2 j . Since the anti- 
bound state is separated by a gap of 0( J) from the two- 
magnon continuum at q 7r, we may use (<P(2,q) | V 1 ^ q)) ~ 



1 in that region, leading to 



\{^l,q)\fr { 2,q))f 



X ^-[%,)-%9)])' 



(14) 



Using Eq. ([10]) and introducing the frequency cu q 
J/2 (1 — cosg) this can be rewritten as 
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Since momentum enters only through u> qi we may intro- 
duce the corresponding density of states and replace the 
sum by an integral with respect to uj q . This results in 



C 



AAA 



j e 0(2w-B*) 



2tt y/J/uj - 1 



e(w)e(j-w), (i6) 



where Q(lj) is the Hcavyside function. Fourier trans- 
forming this to the time domain we get 



il_ e -P(B*-J) e *(B*+J/2)t 



(17) 



where 2o,i(z) are modified Bessel functions of the first 
kind. At low temperatures, i.e. for (3 J 3> 1, we may insert 
their asymptotic forms for \z\ 3> 1, which are Io,i(^) ~ 
e z j\j2n z. Therefore 



Caaa(^) 



i(B* + J)t 

(18) 



y/2(j + lt 

ReCf 



The real part of this, i.e. C AAA (t) = ReC AAA (t), is 
Eq. (4) of our Letter, with lo l = B z . 

To assess the quality of the asymptotic approximation, 
we compare this result with ED for the current autocorre- 
lation function above the saturation field, restricting the 
intermediate-state Hilbert space to two-particle states. 
This corresponds to taking the limit (3 J ^> 1. Due to 
this restriction ED is possible for rather large systems 
sizes N. We choose N = 200. As shown in Fig. the 
agreement is excellent. This validates the approxima- 
tions involved in going from Eq. (fl"0)) to (fl"8|) . 



II. ED Calculations 

In this section we provide supplementary material on 
our ED calculations concerning a potential impact of the 
finite system size. Specifically, we demonstrate that for 
the systems of size N = 18 and time scales tJ < 20, 
as used in the Letter, finite-size effects in the frequency 
shift Su> and the relaxation time r can be neglected. Both 
of these quantities arc determined by assuming a single 
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FIG. 7. The current autocorrelation function C x (t) for (a) < 
t J < 30 and (b) 30 < t J < 60, restricted to the transitions 
into the subspaces M < 1 (one-magnon subspace). Data 
are evaluated numerically by the use of ED (circles) and are 
shown for the parameters f3J — 9, B z / J — 3, and N = 200. 
For comparison, the asymptotic approximation according to 
Eq. (4) of the Letter is included (curves). 



exponential R(t) = exp(— t/r) to model the envelope of 
the decay of the coherent oscillation as in Eq. (3) of the 
Letter. This procedure leads to satisfactory fits to the ED 
results for not too low temperatures, (3J < 3. Examples 
of this are shown in Fig. 2 of the Letter. 

Figure [8] summarizes our results for Suj and r as ob- 
tained from N = 10, 12, 14, 16, and 18 for various mag- 
netic fields at f3J = 2. This figure clearly demonstrates 
that both, Suj and r either display almost no finite size 
dependence or a clear tendency towards saturation as a 
function of 1 /N. In all cases shown, the absolute changes 
in going from N = 16 to 18 are negligible. Significant fi- 
nite size effects can only be seen at systems sizes N < 14. 
Therefore, for the temperatures and fields considered in 
the Letter, it is justified to use TV = 18 data to obtain 
Suj and t (sec Fig. 3 of the Letter). 



III. DMRG Calculations 

In this section we provide details of our DMRG simula- 
tions of the real-time evolution of the transverse current. 
We will focus on three aspects: (i) the coarse graining of 
the current, (ii) the dependence on the initial state and 
its characteristic wave length k, and (iii) the numerical 
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FIG. 8. The extracted (a) frequency shift 5u> and (b) decay 
time r as a function of the inverse size 1/N for the parameters 
/3J = 2 and B z / J = 1, 2, and 3 (symbols). Lines are guides 
to the eyes. 



determination of the frequency shift Suj, depending on 
the maximum simulation time and system size. The sim- 
ulations were carried out with a fixed discarded weight, 
which we have varied from 10 -4 to 10 -6 in order to check 
convergence of our results. 

Figure |9] shows the Fourier spectrum of the transverse 
current for two different coarse graining schemes. The 
solid, black line is taken from Fig. 4(c) of the Letter where 
we average the current over the left half of the chain 
(labeled by a = N/2 in Fig. |9j). The dashed, red line 
is the result from averaging the current over five sites 
in the middle of the chain (denoted by a = 5); in this 
example, the coarse graining is taken over sites r = 30, 
31, 32, 33, 34 in a system with N = 64. As the figure 
clearly shows, the position of the two maxima, i.e., the 
one at the Larmor frequency ujl and the one at ujl + Suj, 
as well as their (relative) weights are insensitive to the 
coarse graining. 

Turning to the initial states of the real-time evolution, 
they are constructed by adding a site-dependent trans- 
verse field via a term J2 r B^S* to the Hamiltonian with 
= B x cos(27rfcr/7V). Figure [TU] illustrates that a cen- 
tral result of our Letter, namely the positions of the max- 
ima in the current's Fourier spectrum, does not depend 
on the wave vector k of the perturbing field B% in the 
long wave-length limit. Computational constraints lead 
to a decrease of the maximum simulation times accessi- 
ble as k is increased. For a comparison of k = 1/2 and 1 
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FIG. 9. Coarse graining: The solid, black line is the 
frequency-dependent current averaged over the left half of the 
system while the dashed, red line is averaged over five sites, 
counting away form the center of the system with TV = 64 (to 
be precise, these are sites r = 30, 31, 32, 33, 34). The main 
features, i.e., the dominant frequencies do not depend on the 
coarse graining. 
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FIG. 11. Numerical determination of lul + Sto in dependence 
of the simulation time, for three different system sizes: N = 
32 (solid, black line), N = 48 (dash-dotted, blue line) and 
N = 64 (dash-dotted, red line). The final values are all within 
the overall numerical accuracy and no systematic finite size 
effects are visible. The inset shows the spectra at tmaxJ = 80 
for N = 32 (solid, black line) and N = 64 (dashed, red line). 



this implies that we have to confine ourselves to t J < 40. 
Therefore, the features in Fig. [TU] are broader than in 
Fig. 4 (c) of our Letter. 

In Fig.[TTJ we address the convergence of the position of 
the second maximum in i^/ 2 ( w ) at lul +Slo as a function 
of the simulation time and the system size. Since we per- 
form a discrete Fourier transform on a finite time window 
to obtain i^ 2 ( a; ) f rom the real-time data i^y 2 (£)) the 
resulting spectrum depends on the maximum simulation 
time. After a minimum time needed to resolve the two 
frequencies has been reached we compute the spectrum 
after each time step and extract both frequencies. The 
quantity displayed in the figure is u>l + Sco averaged over 




FIG. 10. Dependence on the initial state: Spectrum of the 
transverse current for two different values of k: k = 1/2 
(solid, black line) and k = 1 (dashed, red line). In both cases, 
B x /J = 1, B z /J = 1 and N = 64. The Fourier transform is 
taken with t ma _ x J — 40 due to the higher numerical costs at 
larger k. 



these spectra obtained from t ma ^J = 20 up to the max- 
imum simulation time t max . For N = 32 and 48 (solid, 
black and dash-dotted, blue line, respectively), ujl + 5lo 
weakly increases with i max , while for N = 64 (dashed 
red line), the convergence is much faster. The finite-size 
effects in lu^ + Suj are non-monotonous and result in a 
small uncertainty of about 2-3%, well within the overall 
numerical accuracy of the simulations. The inset shows 
the spectrum of the current as a function of frequency at 
t max J = 80 for N = 32 and N = 64. While the dom- 
inant frequencies are only slightly affected by finite-size 
effects, spurious additional peaks appear on the smaller 
system that are irrelevant for the results and discussion 
presented in the Letter. 



IV. Comment on Units 

In the Letter all quantities are expressed in units of 
the exchange coupling constant J. Moreover, and to ab- 
breviate the theoretical analysis, Planck's constant h has 
been set to unity as usual. Similarly, the Bohr magne- 
ton /xb and the spin Lande factor g$ have been absorbed 
into the definition of the magnetic field. The correspon- 
dence between these units and Sl-units is summarized in 
Tab. HI Therein, J is expressed in units of temperature, 
i.e., divided by the Boltzmann constant k-Q. 
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quantity 


unit 


Letter 


SI 


magnetic field 


J 


ks/ifJ-B gs) (J/fe) 
« 0.744 T/K (J/fc B ) 


frequency 


J 


fc B /^ (J/fe) 

« 1.309 • 10 n Hz/K (J/fc B ) 


time 


1/J 


ft/fcs (fea/J) 

« 7.638 • 10- 12 sK (fee/ J) 



TABLE I. Correspondence between the units in the Letter 
and Si-units. 
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